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It has been recently shown that puzzling excess events observed by the LSND and MiniBooNE 
neutrino experiments could be interpreted as a signal from the radiative decay of a heavy sterile 
neutrino (i>h) of the mass from 40 to 80 MeV, with a muonic mixing strength |l7 M h| 2 ~ 10~ 3 — 10~ 2 , 
and the lifetime 10~ n < th < 10 -9 s. We discuss bounds on |i/ M /t| 2 obtained from the recent 
precision measurements with muons and show that they are consistent with the above limits. If 
the Vh exists its admixture in the ordinary muon or kaon decay would result in the decay chain 
H eveVh — s- eue'yv or K — y [Wh M7^i respectively. We propose a new experiment for a sensitive 
search for these processes in fi and K decays at rest allowing to either definitively confirm or exclude 
the existence of the v^. To our knowledge, no experiment has specifically searched for the signature 
of radiative decay of massive neutrinos from particles decays as proposed in this work. The search is 
complementary to the current experimental efforts to clarify the origin of the LSND and MiniBooNE 
anomalies. 

PACS numbers: 14.80.-j, 12.20.Fv, 13.20.Cz 



I. INTRODUCTION 

Over the past 10 years there is a puzzle of the 3.8 
a event excess observed by the LSND Collaboration 
[l|]. This excess originally interpreted as a signal from 
— > v e oscillations was not confirmed by further mea- 
surements from the similar KARMEN experiment 
The MiniBooNE experiment, designed to examine the 
LSND effect, also found no evidence for — > v e oscilla- 
tions. However, an anomalous excess of low energy elec- 
tronlike events in quasielastic neutrino events over the ex- 
pected standard neutrino interactions has been observed 
[3j . Recently, MiniBooNE has reported new results from 
a search for — > T7 e oscillations [I[ . An excess of events 
was observed which has a small probability to be identi- 
fied as the background-only events. The data are found 
to be consistent with — > v e oscillations in the 0.1 eV 2 
range and with the evidence for antineutrino oscillations 
from the LSND. 

In recent work d| (see also @, Q) it has been shown 
that these puzzling results could all be explained in a con- 
sistent way by assuming the existence of a heavy sterile 
neutrinos (yh)- The i>h is created in neutral-current 
interactions and decay subsequently into a photon and a 
lighter neutrino v in the LSND and MiniBooNE detec- 
tors, but it cannot be produced in the KARMEN exper- 
iment due to the high energy threshold. The Vh could be 
Dirac or Majorana type. The Vh could decay dominantly 
into a 71/ pair if, for example, there is a large enough 
transition magnetic moment between the Vh and v mass 
states. Assuming the Vh is produced through mixing with 
the combined analysis of the LSND and MiniBooNc 
excess events suggests that the vh mass, mixing strength, 
and lifetime are, respectively, in the range 

40 < m h < 80 MeV, 10" 3 < \U„ h \ 2 < 10~ 2 , 

10~ n < t h < 10- 9 s. (1) 

A detailed discussion of consistency of these values with 



the constraints from previous searches for heavy neutri- 
nos § is presented in Briefly, the mixing of (P) is 
not constrained by the limits from the the most sensi- 
tive experiments searched for extra peaks in two-body 
7T, K decays ||, because the Vh mass range of |T]) is out- 
side of the kinematical limits for 71^,2 decays, and not 
accessible to K^2 experiments due to experimental res- 
olutions. The parameter space of ([1]) cannot be ruled 
out by the results of high energy neutrino experiments, 
such as NuTcV or CHARM, as they searched for Uh's of 
higher masses (rrih > 200 MeV)_ decaying into muonic fi- 
nal states (^7w, nfiv, /lev, ..) |8j, which are not allowed 
in our case. The best limits on |J7 M ^| 2 derived for the 
mass range (|TJ) from the search for —¥ e + e~v decays 
in the PS191 experiment [9(, as well as the LEP bounds 
[l(J, are found to be compatible with ([TJ assuming the 
dominance of the Uh — > 7^ decay. New limits on mix- 
ing lU^h] 2 obtained by using the recent results on pre- 
cision measurements of the muon Michel parameters by 
the TWIST experiment [ill are a ^ so f° un d to be con- 
sistent with |T]). Finally, the most stringent bounds on 
|C//i^| 2 coming from the primordial nucleosynthesis and 
SN1987A considerations, as well as the limits from the 
atmospheric neutrino experiments, are also evaded due 
to the short Vh lifetime. 

As mentioned above, the most natural way to allow 
the radiative decay of heavy neutrino is to introduce a 
nonzero transition magnetic moment (ut r ) between the 
Vh and v mass states; see e.g. [I2j[I3|- Such coupling 
of neutrinos with photons is a generic consequence of 
the finite neutrino mass. Observations of the neutrino 
magnetic moment could allow us to distinguish whether 
neutrinos are of the Dirac or Majorana type since the 
Dirac neutrinos can only have flavor conserving transi- 
tion magnetic moments while the Majorana neutrinos 
can only have a changing one. In addition, Dirac neu- 
trinos can have diagonal magnetic moments while Majo- 
rana neutrinos cannot. The nonzero magnetic moment of 



the neutrino, although tiny, is predicted even in the stan- 
dard model (SM) . The detailed calculations of the radia- 
tive neutrino decay rate in terms of the neutrino masses 
and mixings were performed long ago, see e.g. EMI. 
The radiative decay mode could even be dominant, if the 
fitr value is large enough; see [H, [l3| . Originally, the 
idea of a large (Dirac) magnetic moment (> 10~ n /is, 
where /is is the Bohr magneton) of the electron neutrino 
has been suggested in order to explain the solar neutrino 
flux variations (l7| . Taking into account that in many 
extensions of the standard model the value of the fi tr is 
typically proportional to the Vh mass, the intention to 
make the radiative decay of a Vh S 100 MeV dominant 
by introducing a large transition magnetic moment (or 
through another mechanism) is not particularly exotic 
from a theoretical viewpoint. Such types of heavy neu- 
trinos are present in many interesting extensions of the 
standard model, such as GUT, superstring inspired mod- 
els, left-right symmetric models and others, for a review; 
see e.g. Ref . |12| . 

The Vh — > "/v decay rate due to a transition moment 
fitr is given by [3 
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where m„(< mj,), is the mass of the lighter neutrino in 
the Vh — > iv decay, /is = ^ s ^ ne Bohr magneton, 
and m e is the mass of electron. Assuming m v -C m^, 
and Th < 10~ 9 , which came from the requirement for the 
Vh — > iv decays to occur mostly inside the MiniBooNE 
fiducial volume, results in Q 



> 3.7 x 10- 



(3) 



A detailed discussion of the interpretation of the Vh — > "/v 
decay in terms of transition magnetic moment and consis- 
tency of the value © with the constraints from previous 
experiments is presented in [5j (see, Sec. VI. F). Let us just 
note, that the v mass state should not necessarily be the 
dominant mass state of an ordinary neutrino, v e , v^ or 
v T . It could be a sterile neutrino as well. In this case, 
there are no any constraints on the transition magnetic 
moment between two sterile states at all. There are also 
no stringent experimental bounds on the fitr between Vh 
and v, if v = v T , because there are no intense r neutrino 
beams. First experimental limits on /i tr between a heavy 
sterile neutrino and muon neutrino have been obtained 
in Ref. [li| , based on the idea of the Primakoff conversion 
v^Z — > VhZ of the muon neutrino into a heavy neutrino 
in the external Coulomb field of a nucleus Z , with the 
subsequent Vh — > jv decay. By using the results from the 
NOMAD experiment [20|], a model- independent bound 
Mtr i5 4.2 x 1CP 8 /U_b was set for the Vh masses around 
50 MeV (see Table 1 and Fig.2 in Ref.[l|), which is also 
consistent with Eq.Q. 

If the Vh is indeed a component of v^s, it would be 
produced by any source of v^ according to the proper 
mixing and phase space and helicity factors [HI, [2l[ . In 




FIG. 1: Schematic illustration of the production and subse- 
quent radiative decay of heavy neutrinos in the ordinary muon 
decay. 



particular, the Vh could be produced in charge-current 
weak interactions of muons or kaons. For example, for the 
mass range of ([1} the nonzero mixing | U^h | 2 would result 
in the decay \x —> ev e Vh- The muon, which is normally 
decays into an ev e pair and a v^, might instead decay 
to an ev e pair and a heavy neutrino Vh which decays 
subsequently into "/v, as schematically illustrated in Fig. 

ni 

In this work we show that the recent precision mea- 
surement results obtained with muons are consistent with 
(P) , and that the puzzle of the LSND-MiniBooNE excess 
events could be uniquely resolved by the proposed new 
experiment on searching for the muon and/or kaon de- 
cay chains, fi — > ev e Vh — > ev e jv and K — > \n>h — > ^"fv, 
respecively. 



II. CONSTRAINTS FROM MUON PROCESSES 

For completeness of the analysis of constraints reported 
in [5J], let us examine first the recent precision measure- 
ments results obtained with muons. Note that for the 
mixing and mass regions of |T]) , the inclusion of the heavy 
neutrino effect results in the /i — > ev e Vh decay rate which 
can be well approximated by [2l|, [22[ 
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and in the corresponding branching fraction 

B(fi ev e v h ) ~ 10~ 5 - 10~ 3 , (5) 

which is in the experimentally accessible range. More de- 
tailed calculations of the — ¥ ev e Vh decay rate including 
arbitrary Vh weak couplings can be found in [23| . 



3 



Muon lifetime 



in the 2a limit B((i -> ev e v h ) < (2.6-4.1) x 10" 4 , which 
is consistent with (JSJ) - Interestingly, for the large values 



Very recently the MuLan Collaboration has reported of which are expected for the fi 



ev e v h 



eveiv 



on measurements of the mean lifetime r M of positive 
muons to a precision of 0.6 ppm [24j . Using the new world 
average, r M = 2.1969803(22) fis and the relation between 
the muon lifetime and the Fermi constant Gp t^ 1 = 

192^3* (1 + wnere A is the sum of phase space, QED 
and hadronic corrections, results in new determination of 
the Fermi constant G F = 1.1663788(7) x 10~ 5 GeV~ 2 to 
a precision of 0.6 ppm 24]. The mixing of the vy L into 
Vfi would decrease the determined value of G F . To esti- 
mate the allowed contribution from the /i — > ev e Vh decay, 
one could compare the experimentally measured muon 
decay rate to a predicted one, by using G' F extracted 
from another measurements which are not directly af- 
fected by the contribution from heavy neutrino. One 
possible way is to use the pure leptonic decay rate of 
the tau r(r — » ev e v T ), which provides the correspond- 
ing Fermi constant G T F e = 1.1668(28) x 10~ 5 GeV" 2 
Q. Comparing it with MuLan values for Gp one finds 
AG F = GV-G*i, < 5 x 10- 3 GeV- 2 (90%C.L.). Taking 
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into account (J4j) leads to the bound \U)j,h\ 2 < 8 x 10 
which is consistent with (p}. One can also use a num- 
ber of indire ct p rescriptions for extracting of precise val- 
For example, one can define G' F — 
where @w> rriz and Ar are the 



ues of Gf [25 



4-Ka 



sm 2 20w (mz)(l — Ar) 

Weinberg angle, the mass of the Z gauge bosons ex- 
tracted from the precision measurements at LEP, and 
a factor for radiative corrections, respectively. Using the 
values of Ow> m z an d Ar reported in (25|, one can ob- 

Comparing it with G F and adding statistical and system- 
atic errors in quadrature, one finds at 90% CL, AGf — 
G' F -G F < 4.1 x 10- 3 GeV- 2 (90%C.i.), which leads 
to the bound |t^h| 2 < 7 x 10~ 3 , which is also consistent 
with ([1]). 



tain G' F = 1.1672(±0.0008) 



B. Rare muon decays 

As the final states of the decay /i — > ev e Vh — > ev e ^v 
and the radiative muon decay /x — > ev^^ are iden- 
tical, the decay rate of the former can be constrained 
from the precise measurements of the branching fraction 
B(fi — ¥ eVeV^) of the later. The most precise measure- 
ments of the radiative muon decay are reported by the PI- 
BETA Collaboration [26[ . The measured branching frac- 
tion B eX p(/j, — > eVeV^) = (4.4 ± 0.1) x 10 -3 is in a good 
agreement with the predicted value Bsm((J< — > eVeV^) — 
4.3 x 10~ 3 for the photon energy E 1 > 15 MeV, and the 
e — 7 opening angle O e7 > 45°. Thus, the contribution 
from the decay /i — > e.v e Vh is allowed to be at the level 
\B exp (ii ev e v^) - B SM {p. -> ev e v^)\ ~ (1.0 ± 1.0) x 
10 -4 , which, taking into account the efficiency for the 
above selection criteria, results for the Vh mass range ([T]) 



events, PIBETA has observed an excess of misidentified 
nonradiative events attributed to the so-called "splash- 
back" events [26|]. The size of this effect is comparable 
with . It was assumed that these rare events originated 
probably from Michel decays when the shower compo- 
nent of the decay electron entering an electromagnetic 
calorimeter is emitted at a very large angle with respect 
to the primary electron momentum direction. 

The decay chain /i — > ev e Uh — > ev e ^v could also con- 
tribute to the background for the lepton flavor violating 
decay /i — > ej whose branching fraction is constrained to 
be B exp (fi -> e-f) < 1.2x 10 -11 by the MEGA experiment 
[2^ | . In this experiment, to avoid background from the 
radiative muon decay, only back-to-back ej pairs were 
selected for analysis. The energy of the electron was re- 
quired to be around the end point of E e = 52.8 MeV 
within the energy resolution of 0.23 MeV [27[. Because 
the maximal allowed electron energy in the decay fi — > 



ev e Vh is significantly less, 



m. 



MeV, this potential background decay mode was rejected 
by the above selection criteria. 



C. Muon capture 

Recently, bounds on lE/^fcJ 2 have been obtained from 
the measurements of the radiative muon capture rate on 
hydrogen /xp — > v^wy, for the 7 energy threshold of E 1 > 
60 MeV [28| (see also 29]). The new limits seem are 
in tension with |C/ M /i| 2 values from (JTJ) , although their 
exclusion strength is a factor of few. For Uh masses from 
40 to 70 MeV the limits are, respectively, from \U a h\ 2 < 
7 • 10- 4 to IIVJ 2 < 5 • 10- 4 (Fig. 4, a = 1 in @), 
while the corresponding 2a lower bounds of the allowed 
LSND-MiniBooNE parameter region (Fig. 24, a = — 1 in 
U) are in the range from lU^] 2 > 2 ■ 10" 3 to \U^ h \ 2 > 
3T0~ 3 (InRef. 28] the definition of sign of the asymmetry 
parameter a of the photon angular distribution in the rest 
frame is opposite to the one in H). However, the more 
conservative 3 a LSND-MiniBooNE lower bounds for the 
same mass range are calculated to be, respectively, from 
\U^ h \ 2 > 4 ■ 10~ 4 to \U^ h \ 2 > 5 ■ 10~ 4 , and the tension 
is not apparent. Furthermore , regardless of the level 
of the exclusion strength, there is also a direct way to 
evade the radiative muon capture limit by shifting the 
photon energy spectrum from the i>h — > decay towards 
the lower energy region below 60 MeV by allowing the 
neutrino v to be massive enough, as suggested in [28j |. 



III. EXPERIMENTS TO SEARCH FOR THE 
u h -> 7^ DECAY 

In order to (dis)prove the interpretation of the 
LSND-MiniBooNE excess events, we propose a new ex- 
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FIG. 2: Schematic illustration of the proposed beam-dump 
type experiment to search for the decay fi + — > e + v e Vh — ► 
ev^v. Shown are the beam defining scintillator counter 
S, the active target T surrounded by the electromagnetic 
calorimeters ECAL1 and ECAL2, and the ECAL1 endcap 
counter EC used as a light guide for the light produced in 
the target. The source of muons is the n — > /.t + v decay at 
rest in T. Electrons and photons from the ordinary muon de- 
cays n + — > e + 1^1^(7) are dumped in the assembly, consisting 
of the ECAL1 and a lead absorber Al, used as a shield for 
the ECAL2. Electrons and 7's from the n, /i decays in flight 
before the S counter are absorbed by the shield A2. The Vh 
produced through the mixing with the muon neutrino pen- 
etrates the (ECAL1+ disassembly and decays in flight in 
the free space between the ECAL1 and ECAL2 into a light 
neutrino and photon, which is detected by the ECAL2. The 
insert shows the front view of the ECAL1 area. 



periment to search for the Vh — > jv decay chain from the 
muon or kaon decays at rest with a sensitivity in |J7^^| 2 
several orders of magnitude higher than the |t/ M /i| 2 values 
of®. 



A. Search for the jj, — > ev^vu — > eVefyv decay 

The main components of the detector are schemati- 
cally shown in Fig. [5J A beam of positive pions with 
momentum of ~ 70 MeV/c, e.g. from the 7rM3 pion 
beam line at PSI, , defined by a scintillator counter S, 
is stopped in an active target (T) instrumented with en- 
ergy deposition and time readout. The source of muons 
is the 7r — > u + v decay at rest in T. The target T is 
surrounded by an almost 4tt hermetic shield consisting 
of an electromagnetic calorimeter ECAL1 and a lead ab- 
sorber (Al), where all the electrons and photons from 
the tt,u decays in T are dumped. The ECAL1 is sur- 
rounded by another ECAL2 to detect photons from the 
Vh — > "fv decays in flight. It is assumed that the Vh pro- 
duced via the mixing in the fj, — > ev e Vh decay is a weakly 



interacting particle, which penetrates the (ECALl+^41) 
assembly without significant attenuation and decays sub- 
sequently in the free space between the calorimeters, as 
shown in Fig. [2] Electrons and 7's from the 7r, u decays 
in flight before the S counter, cannot hit the ECAL2, as 
they are absorbed by the shield A2. The experimental 
signature of the decay chain /1 — > ev e Vh — > ev e jv is two 
signals of the energy deposition in the calorimeters sep- 
arated in time by an interval At < 1 ns corresponding 
to the time-of- flight. The physical background for the 
signal events is expected to be very small, because only 
neutrinos can penetrate the shield, but they produce a 
negligible number of events in the ECAL2. As the pro- 
posed experiment is of the beam-dump type, there is no 
special requirements for the purity of the beam. The 
amount of muons and electrons from n decays in flight 
could be comparable to the number of pions at the end 
of the beam line. 

To estimate the sensitivity of the proposed experiment 
a feasibility study based on simplified Monte Carlo simu- 
lations, similar to the one described in [30| . has been per- 
formed. The pions are stopped in the T, which is a plastic 
scintillator with a diameter of 5-10 mm and a height of 
10 mm. According to simulations the decay muon came 
to rest passing about 1 mm in T and depositing about 
4.2 MeV in the T, which can be used for its identifica- 
tion. The ECAL1, is an array of 7 BGO counters, as 
schematically shown in Fig. [2j each of 55 mm in diame- 
ter and 220 mm long, which were previously used in the 
PSI experiment on precise measurements of the 7r e2 decay 
rate [31j. The ECAL2 is an array of the same counters 
with electron energy resolution of ~ 2.5% at 50 MeV. 
The readout of the energy deposition in the ECAL2 is 
triggered by a tag signal of the electron appearance from 
the decay chain tt — > a — > e + anything, which is defined 
by a coincidence of a signal from a stopped pion, a de- 
layed signal from the stopped decay muon and a delayed 
signal from the EC ALT The light signals produced in 
T could be readout through the ECAL1 endcap crystal 
(EC) which acts as a light guide as shown in Fig. [2] The 
T signals could be distinguished from the EC signals due 
to their significantly different decay times by using the 
technique described in detail in Ref.|32]. 
The significance of the Vh discovery in the proposed 
experiment scales as S = 2(^/n s + rib — y/rib), where 
n s , njj are the number of detected signal and background 
events [33|. The number of events expected from the 
fi — > ev e Vh — > ev e -^v decay chain is calculated as 

n s ~ n^B^i -> ev e v h )B(v h -> -fv)P Uh ^ lv f 1 e 1 t, (6) 

where n^ is the muon stop rate in the target, B(vh — > 
1V ) ~ 1 i, Pvu^p — exp(—l I 'cjTh) is the probability 
for the Vh to decay in flight in free space , where I is an 
effective (ECAL1+A) thickness and 7 ~ 1 is the gamma- 
factor of the Vh, fj > 0.8 is the fraction of events with 
the energy deposition in the ECAL2 > 10 MeV, e 7 ~ 0.1 
is the average 7 detection efficiency, and t is the running 
time. Simulations show that the energy spectrum of de- 
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cay photons is well above ~ 10 MeV for the Vh masses 
of JT]) and a wide range of the v masses. To estimate 
rib, the following background sources are considered: (i) 
the radiative 7r//x decays, which have the branching frac- 



tion 



1(T 2 for £L > 10 MeV 



To suppress this 



background the effective thickness of the (ECAL1+A1)- 
assembly is selected to be I > 16 X$ resulting in < 10~ 8 
decay photons per pion stop penetrating the shield with- 
out interactions; (ii) bremsstrahlung of decay electrons in 
the ECAL1 results in a leak of low energy photons to the 
ECAL2, which can be rejected by the energy cut E 7 > 10 
MeV; (hi) accidental 7's from 71", fi radiative decays in 
flight are rejected by the shield A2 and by the require- 
ment to have only one entering particle for a particular 
event. The number of 7's from the [i — > ev e v^ decays 
scattered in the S counter is found to be small; (iv) ac- 
cidental coincidences from cosmic rays can be removed 
by an active veto and can be neglected; (v) the neutron 
background is expected to be small, as the setup could 
be located at a large (> 30 m) distance from the proton 
target. Finally, we found that the expected background 
rate rib i$ 10 _8 /7r decay is dominated by (i), and is well 
controlled by the choice of the absorber Al thickness. 
Assuming S > 5 and the muon stop rate n M ~ 3 • 10 4 p/s 
one could expect the sensitivity in the 11 — » ev e Vh decay 
branching ratio as small as B{p — > ev e Vh) < 10~ 9 for the 
beam exposure ~ 1 month. 

In the case of signal observation, the value of the Vh 
mass can be evaluated from the end point of the 7-energy 
spectrum. To cross check the origin of signal events, one 
could use an additional shield in the space between the 
ECALs to absorb a fraction of photons from the Vh — > iv 
decays in flight, which can be well predicted. If no signal 
events are observed, the limit on the mixing strength: 



\U»h\ 2 < 10- 8 e°- 3 /^M 



m,, 



-7/2 



(7) 



could be set as a function of the Vh mass and lifetime. 
Using (JTJ) , one can see that for the mass range of (U} and 
lifetime values in the range 5 • 10~ u < 77, < 10 -9 s the 
limits on the mixing strength are |J7 M ft,| 2 < 10~ 7 — 10~ 4 
and hence, the |i/ M /i| 2 values of |T|) would be firmly ex- 
cluded. The choice of the (ECALl+Al)-assembly thick- 
ness compromises the rejection factor of background 7's 
from the [i — > ev e v^ decay and the sensitivity in |J7 M ^| 2 . 
For the lifetimes as short asT^<5-10 _11 s the vast 
majority of v^s decays in the vicinity of the target and 
the limit is less restrictive. To improve sensitivity for 
this lifetime region, one could replace the ECAL1 with 
an assembly of thin plastic scintillator counters to detect 
decay electrons and the Al with an absorber made of a 
higher- Z material, e.g. tungsten, such that the overall 
thickness of the shield is reduced to a few centimeters. 



B. Search for the K — > (j,Uh — ► \IT(V decay 

A similar search for the i>h can also be performed 
with another experiment, analogous to the one discussed 
above. In this experiment we propose to search for 
the K — > ixvh —> l^iv decay chain from kaon decays 
at rest. The main component of the detector are the 
same as shown in FigJ5J but the incident beam is com- 
posed mainly of positive kaons. The tag for the two-body 
K — > (iv decay is the energy deposition around maximal 
value of 150 MeV in the ECAL1 from the stopped muon 
kinetic energy, which is expected for the neutrino masses 
from to 80 MeV. Similar to the above muon experiment, 
the experimental signature of the K — > (iv^ — > fi^v decay 
is two signals of the energy deposition in the calorime- 
ters ECAL1 and ECAL2 separated in time by the time 
interval At < 1 ns corresponding to the Vh time-of-flight. 
The number of K — > \Wh events is defined by the mix- 
ing lUfih] 2 and by the phase space and helicity factors 
which depend on the mass [3J]. For the mass interval 
m Vh ~ 40 — 80 MeV the chirality-flip is mostly due to the 
sterile neutrino mass which results in 

2 

T(K -> iiv h ) « T{K -> ^)\U^p(-^) , (8) 



where term py -^r) includes both phase space and he- 
licity factors [3J]. Using ([8]). we find that the branching 
fraction of K — > i±V} x is in the experimentally accessible 



range: 



Br(K -> fiv h ) ~ 10" 3 - 10" 



(9) 



for heavy neutrino masses in the range 40 - 80 MeV. As- 
suming the kaon stopping rate uk — 3 TO 4 fi/s one could 
expect the sensitivity in the K — > fxi/h decay branching 
ratio as small as B(K — > fii/h) S 10~ 9 for the beam expo- 
sure ~ 1 month. The expected sensitivity of the kaon ex- 
periment is comparable with ([7)1, however for the Vh life- 
time values < 10 -10 s it is even higher due to the higher 
energy of iVs, and, thus the longer vh decay length. This 
gives an advantage for searching for the short-lived v^ L in 
K decays compare to the muon experiment. 



C. Search for the fi — s> ev e vu ~> ev a ^v decay with 
cosmic muons 

Finally, note that a search for the /i — > ev e Vh — > eu e ji/ 
decay events may also be possible with existing cosmic 
muon data accumulated by large Cherenkov neutrino de- 
tectors, such as SuperK MiniBooNE [II], etc.. The 
idea is to search for an excess of e — 7 events from stopped 
cosmic muon decays, which are characterized by a clean 
Cherenkov ring from the electron from muon decay and 
the well separated ring from a Compton-scattered elec- 
tron from the gamma from heavy neutrino decay [3^ |. 
The photon spectrum from the radiative muon decay 
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drops quckly with photon energy (~ 1/E^), therefore, 
one would expect a low background level in the signal 
region of E 1 > 20 — 30 MeV and large opening angle 
O e7 > 30°. However, as the i>h decay length (< 30 cm) 
is significantly less than the effective photon Compton- 
scattering length ( ~ 60 cm) , both the signal and " splash- 
back" background events, discussed in Sec. II B, are ex- 
pected to have similar spatial distributions with respect 
to the primary decay vertex. Thus, one cannot effec- 
tively suppress these background events by using spatial 
cuts. The use of cuts only on the photon or the sum of 
the electron and associated photon energy deposition can 
not be efficient enough, in particular, due to a moderate 
energy resolution of Cherenkov detectors in the 10 - 50 
MeV energy range, e.g. ~ 15% at 53 MeV in MiniBooNE 
[36| . An estimate of the sensitivity of such an experi- 
ment requires detailed simulations, which itself might be 
inaccurate in this extreme "splashback" regime. This is 
beyond the scope of this work. 

A search for fi — » ev e Uh — » ev e piv events could also 
be performed by the ICARUS [3g. In this experiment 
the decay electrons energy spectrum from a sample of 
stopping cosmic muon events recorded during the test 
run of the ICARUS T600 detector was studied. The 
detector allows the spatial reconstruction of the events 
with fine granularity, hence, the precise measurement of 
the range and dE/dx of particles with high sampling 
rate. The measured energy (E exp ) resolution for elec- 
trons below 50 MeV is finally extracted from the sim- 
ulated (Egim) sample, obtaining (E exp - E sim )/E sim = 
2% + ll%/E°- 5 [MeV\. The detector has good spatial 
and energy resolution, and is capable of distinguishing 
between electron and photons, hence, a better sensitiv- 
ity for the search for \i — > ev e Vh ~^ zveiv events can be 
expected from the analysis of recorded data. 

Very recently, another interesting idea of searching for 
jVs has been proposed [39]. It has been noticed that a 
copious production of i/^'s from K decays in flight could 



occur in air showers. Taking this into account, a sensitive 
search for signals from subsequent Vh radiative decays in 
neutrino telescopes, like ANTARES or the DeepCore in 
IceCube, could be performed. 



IV. SUMMARY 

In summary, we showed that the recently proposed ex- 
planation of the puzzle from the LSND, KARMEN and 
MiniBooNE experiments in terms of the radiative decay 
of a 40-80 MeV sterile neutrino, could be uniquely probed 
by the proposed new experiment on direct search for this 
decay in the muon and/or kaon decays at rest with the 
sensitivity in a corresponding branching fraction as small 
as a few parts in 10 9 . The quoted sensitivity could be ob- 
tained with the proposed setup optimized for several of 
its properties, such as the absorber Al thickness, during a 
month of data taking and would allow us to either confirm 
or rule out the existence of the Vh- Two different designs 
could be implemented for muon experiment to cover the 
range of the lifetime values 10 -11 — 10~ 9 . The proposed 
searches are complementary to the current experimen- 
tal efforts to clarify the origin of the excess events ob- 
served by the LSND and MiniBooNE experiments. This 
enhance motivation for the proposed experiments to be 
performed in the near future. To our knowledge, no ex- 
periment has specifically searched for the signature of 
radiative decay of massive neutrinos from muon or kaon 
decays as proposed in this work. The search for events 
from the Vh — ► iv decay can also be performed by using 
the existing cosmic muon data, collected b y n eutrino de- 
tectors such as SuperK [H|, MiniBooNE [36j|, ICARUS 
[11], as well as in neutrino telescopes (39|. 
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